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We present a comprehensive theoretical study of thermodynamic properties of superconductors 
with a dilute concentration of magnetic impurities, with focus on how the properties of the supercon¬ 
ducting host change if the magnetic moments of the impurities order ferromagnetically. Scattering 
off the magnetic impurities leads to the formation of a band of Yu-Shiba-Rusinov states within the 
superconducting energy gap that drastically influences superconductivity. In the magnetically or¬ 
dered system, the magnetization displays a sudden drop as function of impurity density or magnetic 
moment amplitude. The drop occurs as the spin-polarized impurity band crosses the Fermi level 
and is associated with a quantum phase transition first put forward by Sakurai for the single impu¬ 
rity case. Taking into account that the background magnetic field created by the ordered impurity 
moments enters as a Zeeman shift, we find that the superconducting phase transition changes from 
second order to first order for high enough impurity concentration. 

PACS numbers: 74.78.-w, 74.78.Na, 74.20.Fg, 74.62.Dh, 73.63.-b 


I. INTRODUCTION 

In conventional and most known high-temperature su¬ 
perconductors Cooper pairs are formed of electrons with 
anti-parallel spins (spin-singlets). On the other hand, 
existence of a ferromagnetic order favors parallel spin 
alignment. Therefore superconductivity and ferromag¬ 
netism are two competing orders and they are known to 
coexist in only a few bulk rare-earth material^® and 
more recently at an interface between two bulk insula¬ 
tors, LaAlOa (LAO) and SrTiOa (STO)P. Though such 
systems are rare in Nature, they can be engineered ar¬ 
tificially in superconducting heterostructures. An exam¬ 
ple of such systems are ferromagnet-superconductor (FS) 
interfaces and SFS Josephson junctions, where a long- 
range proximity effect was observecPHU. Another type of 
system which are a focus of active research are supercon¬ 
ductors with magnetic impurities. The term “magnetic 
impurity” in this context may refer to a single transition 
or rear-earth atorcPtlll or small ferromagnetic islandJ^^. 
Another reason for the current interest in these systems 
is the search for experimental evidence of elusive Majo- 
rana particles. According to theoretical prediction^i^lfi^ 
they emerge at the ends of linear chains of magnetic im¬ 
purities, and are believed to have been observed in recent 
experiment 

From a theoretical point of view, investigation of prop¬ 
erties of superconductors with magnetic impurities dates 
back to the work by Abrikosov and Gorko\E^. Within 
the first order Born approximation (weak impurity scat¬ 
tering) they demonstrated that magnetic impurities lead 
to pair-breaking and within certain parameter regimes of 
the model, gapless superconductivity can emerge. Tak¬ 
ing into account the scattering off a single impurity ex¬ 
actly, to all orders in perturbation theory, later it was 
showrPsHU] that magnetic impurities in a superconductor 
are able to host single-particle bound state^^MUl^ com¬ 


monly known now as Yu-Shiba-Rusinov (YSR) states. 
Each YSR state is spin-polarized and does not have a 
Kramers partner due to explicitly broken time-reversal 
symmetry. As a consequence, occupying or emptying 
such a state, the system undergoes a quantum phase tran¬ 
sition and its ground state changes its parity, at the same 
time gaining or loosing one single-particle spin. Since this 
idea was first put forward by SakuraP^, to our knowledge 
it received a very limited attention resulting in a few 
published theoretical workP^l®^. In these studies they 
usually performed self-consistent tight-binding calcula¬ 
tions with one, two or three impurities to demonstrate 
the phase transition and spin-polarization of the ground 
state. We have to note that YSR states as discussed 
above were found by treating magnetic impurity spins 
as classical. They emerge in a superconductor as an at¬ 
tempt to screen the impurity magnetic moment, similarly 
to the Kondo effect in normal metalJ^. On the other 
hand, treating impurity spins quantum mechanically, it 
is possible to study the interplay of Kondo screening and 
superconductivitji22MSI^ however theoretical treatment of 
such systems becomes much more involved. 

Consider a superconductor with a finite concentration 
of YSR magnetic impurities. An open question is how 
superconducting properties are altered when the impu¬ 
rity spins become ferromagnetically ordered. This can 
happen, for example, due to an external magnetic field, 
impurity exchange interaction (direct or mediated by the 
itinerant electrons) or an intrinsic magnetic anisotropy. 
This issue was studied theoreticalljESHm treating impu¬ 
rity scattering in the Born limit. In this paper we extend 
this study by describing the single-i mpurity scattering 
exactly, using the t-matrix formulatiorP^lRRSl jt is im¬ 
portant to note that previous studies based on the first 
order Born limit took into account only local exchange 
scattering of itinerant electrons off magnetic impurities. 
We argue that, when generalizing the model to the t- 
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(a) 


(b) 



FIG. 1. We consider a superconducting thin film (S) deposited 
on a substrate. Magnetic impurities are homogeneously dis¬ 
tributed within the sample. Magnetic moments of the impu¬ 
rities are either completely unpolarized (a), or ferromagneti- 
cally ordered (b) in the plane of the film. 


matrix approximation, it is necessary to also take into 
account a background magnetic field created by the im¬ 
purity spins in order to obtain physically sensible results. 
This circumstance forbids considering the unitary limit 
of impurity scattering, since increasing the magnitude of 
magnetic moments (unitary limit corresponds to letting 
the moments be infinitely large) drives the system to the 
normal state quickly. In the framework of the qusiclas- 
sical Green’s functions formalisrrPnHSU .^g study the be¬ 
havior of the order parameter and the superconducting 
transition temperature for the case of unpolarized and 
ferromagnetically ordered impurities. For the ordered 
case we find that for high impurity concentrations the 
order of phase transition can change from second to first. 
Furthermore, by computing the magnetization in the sys¬ 
tem (in the ordered case) we are able to demonstrate the 
signature of the quantum phase transition put forward 
by SakuraP^ for a finite concentration of impurities. 

The paper is organized as follows. In Sec. [iT] we 
describe in detail the theoretical model of the system 
and briefly introduce the quasiclassical Green’s function 
formalism. We provide a set of self-consistency equa¬ 
tions and expressions for thermodynamic quantities and 
observables for the two cases: (i) unpolarized and (ii) 
ferromagnetically ordered impurities. In Sec. |III| we 
use the results of the previous section to compute self- 
consistently the transition temperature, order parameter. 


density of states, and magnetization. Section IV contains 
a discussion about the validity of our model and con¬ 
nection to its potential experimental realization. This 
section also contains a summary of our findings and con¬ 
cludes the paper. A few technical details of the calcula¬ 
tion have been collected in the Appendix. 


ties are treated as classical spins, similarly to the model 
proposed by Yu, Shiba and Rusinov^SlEl] ’^g consider 
two cases: impurity spins are (i) randomly oriented or 
(ii) ferromagnetically ordered in the plane of the super¬ 
conducting film. The film thickness is smaller than the 
London penetration depth, so that the magnetic field cre¬ 
ated in the ferromagnetically ordered case does not lead 
to screening currents. The Zeeman shifts can however be 
substantial and we take these into account. We describe 
the impurities by the following Hamiltonians 


N 


= '^i'^0 + avsnij ■ (T)S{r - Vj), 
i=i 


( 1 ) 


N 


= Pnvsm ■ a -t- ^(uq + ausm • ct)5{y - r^), (2) 
i=i 

corresponding to the two cases mentioned above. Local 
scattering off a given impurity at position Vj consists of 
a scalar part parametrized by vq and an exchange part. 
The local exchange scattering is parametrized by a di¬ 
mensionless parameter a, the tunneling amplitude of a 
quasiparticle onto the impurity site, and vs the param¬ 
eter proportional to the impurity magnetic moment!^. 
The unit vector m.j points in the direction of the im¬ 
purity magnetic moment, and cr = (cTa;, (Ty, is the 
vector of Pauli matrices in spin space. For the case of 
unpolarized magnetic impurities, we have “ 0- 

The parameter laj < 1 can take both positive and neg¬ 
ative values depending on the microscopic nature of im¬ 
purities, see Fig. For the case of impurities made 
of transition metal atoms with partially filled electronic 
d-shells, e.g. manganese Mn, the local exchange inter¬ 
action with itinerant electrons in the superconductor is 
anti-ferromagnetic, a > On the other hand, if im¬ 
purities are made of rear earth elements, e.g. samarium 
Sm, with partially filled f-shells, the local exchange scat¬ 
tering is ferromagnetic, a < 0. Finally, for the case of 
ferromagnetically aligned magnetic impurities, there is 
an extra term in the Hamiltonian, see Eq. ([^ . The first 
term in Eq. (§ describes a homogeneous background 
magnetic field created by impurity magnetic moments, 
see Fig. Sb). Here /3 ~ I is a dimensionless fitting pa¬ 
rameter, which has a meaning of a geometrical structure 
factor determined by the actual impurity distribution in 
space, and n is the density of impurities. 


II. THEORETICAL MODEL 

A. Description of the setup 

The system we have in mind consists of an s-wave spin- 
singlet superconducting film with randomly distributed 
magnetic impurities, see Fig. Q. The magnetic impuri- 


B. Quasiclassical Green’s function 

In our calculations we use the quasiclassical Green’s 
function formalisnP2HIII_ Since we aim at describing equi¬ 
librium properties of our system, the central object of 
interest is the Matsubara Green’s functioiP^ 5 (e„, Pf, r)- 
The “hat” denotes a 4 x 4 matrix structure in com¬ 
bined spin and Nambu (or particle-hole) space. The 
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FIG. 2. (Color online) Example of a transition metal ele¬ 
ment (Mn) and a rear earth element (Sm) with their valence 
electron configurations. When placed in a superconductor, 
these elements strip the outer 4s^ and 6s^ shells and behave 
as localized magnetic moments. The Hund’s rule determines 
the type of local exchange interaction (ferromagnetic or anti¬ 
ferromagnetic), which enters the theory through the sign of 
the parameter a. Itinerant electron is depicted by the dotted 
red arrow. 


Eq. (§ as a system of (Riccati-type) transport equations 
for coherence functions, 

[ihwF ■ V + 2 fe „)7 = 7 A 7 + E 7 — 7 E — A, 

, (7) 

{ihvp ■ V — 2 fe „)7 = 7 A 7 + E 7 — 7 E — A. 

The coherence functions for negative and positive Mat- 
subara frequencies are related via a symmetry relation 

7 (e„ < 0 )|e„^-e„ = [ 7 (e„ > 0 )]^ ( 8 ) 

which allows us to express all physical observables in 
terms of only positive Matsubara frequencies. Therefore 
we consider 7 and 7 only for e„ > 0 below. 

2. Analytic continuation: retarded and advanced 
propagators 


Green’s function depends on the Matsubara frequency 
e„ = 2TTkBT(n + 1 / 2 ), the quasiparticle momentum on 
the Fermi surface pf, and the spatical coordinate r. Here 
T is the temperature and ks is the Boltzmann constant. 
The pro pag ator g satisfies the quasiclassical Eilenberger 
equation^ 


[ien% -h,g]+ ih-vp • Vg = 6, (3) 

where fa is the third Pauli matrix in Nambu space, vp is 
the Fermi velocity, and h is the Planck constant. Equa¬ 
tion (§ has to be supplemented by a normalization con¬ 
dition g^ = —TT^i. In Eq. (j^ the self-energy matrix k is 
parametrized as 


Most properties of a system in thermal equilibrium can 
be calculated directly from the Matsubara propagator 
g. However, energy-resolved quantities like, for exam¬ 
ple, density of states need to be computed on the real 
energy axis. So, we need a recipe how to perform an an¬ 
alytic continuation from Matsubara frequencies to real 
energies. In most common cases analytical continuation 
is a formal procedure, which establishes a one-to-one cor¬ 
respondence between the Matsubara g and t he reta rded 
g^ and advanced g^ propagators, which readd^^MI 


9 (^) — > 0)|j£^_).£_|_iO+) 

9 (^) ^ 0) ^e— 20 +■ 


(9) 


For real energies the “tilde”-operation is defined as 


h = 


S A 
A E 


(4) 


y(e,pF,r) = y(-e,-pF,r)*. (10) 


where each element has a 2 x 2 structure in spin space, 
and we have introduced the “tilde”-operation defined as 

y(e„,pF,r) = y(e„,-pF,r)*. (5) 


1. Riccati parametrization 

In order to solve Eq. ( |^ it is convenient to employ 
the Riccati parametrizatioiP^IllIIIl for the propagator g. 
It is realized in terms of “coherence functions” 7 and 7 

^ j53l58l 

» = +=’"(// A)(6) 

g = {l- 77 )-\ F = 07 , 

where T and ± correspond to positive and negative Mat¬ 
subara frequencies, respectively. This allows us to rewrite 


C. Self-consistency equations and observables 

Since the formalism we use is a mean field theory we 
need to supplement the above equations with the corre¬ 
sponding self-consistency equations for the self-energies. 
Because we have no control over the distribution of im¬ 
purities, we have to make several standard yet important 
approximations. In what follows we describe properties 
of the system assuming averaging over different impu¬ 
rity configuration^^. Furthermore we assume a dilute 
concentration of impurities which allows us to use the 
so-called non-crossing t-matrix approximation for find¬ 
ing the impurity self-energj®. The last important as¬ 
sumption, which is usually implicitly used, is that the 
superconducting order parameter Ag (assumed real in 
this paper) is a self-averaging quantity which means that 
in presence of a dilute impurity distribution the order pa¬ 
rameter is homogeneous on the superconducting coher¬ 
ence length scaleS21. After impurity averaging the system 
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becomes homogeneous and thus Eq. Q transforms into 
a system of algebraic equations. 

We introduce for the rest of the paper a set of param¬ 
eters characterizing the impurity subsystem, 

TL 

Uo,S = T^NfVo,S, r = ——, (11) 

TTJV F 

where Nf is the density of states per spin and unit vol¬ 
ume at the Fermi energy in the normal state. Below we 
summarize all the self-consistency equations as well as 
observables considered in this paper for the two models 
of magnetic impurities mentioned above. 


1. Randomly oriented impurities 


We have two types of self-consistency equations: (i) for 
the impurity self-energy and (ii) for the order parameter 
Aq. For the case of unpolarized impurities, see Fig. [^a), 
the coherence functions as well as the superconducting 
pairing self-energy have spin-singlet structure 

7 = *0-270 , 7 = *0-270, As = i(T2Ao. (12) 


Then the impurity self-energy in the t-matrix approxi¬ 
mation can be written as (see Appendix [A| 


^(cra) — |'ao(l ~i~ Uq — a^v^) — i [mq -I- 
+{ul-a^ulf] (25-l)|, 


^imp(^n) — imp (^n) 5 

Ao,imp(e«) = - ^ ho - + (**0 - ^ 70 , 

d = (1 -I- Mq — a^Ug)^ -I- 4a^Mg(20 — 1)^. 


The corresponding “tilded” counterparts Simp and Aimp 
are obtained by applying the tilde-operation, see Eq. ([^. 
It is convenient to introduce two new quantities 


-y. _ -y. 

TT _I _■ ^imp Zjiiup _ A 

^ ' 5 — ^0 


^0,imp- (14) 


It can be shown that in terms of the new quantities the 
solution to the (homogeneous) Riccati equations reads 


7o(en) = i 


Dr, 


Er, + ^Dl + El' 


70 (en) = -70 (Cn), (15) 


and the symmetry relations = E^ and hold. 

Using equations (13l-(15) we can write down a pair of 
self-consistency equations for and 

[ul + a^ul + {ul - a^ul)^] Er^^Dl + El 


Dji — Cn 


Vr, 


Dn — An + r 


[ul - a^ul + {ul - a^ulf]Dny^Dl +El 


Vr, 


Vr, = {l + ul- a^uiY{Di + Ei) + ^a^uiEi. 


(16) 


Let us now find the self-consistency equation for the order 
parameter. By definition, 

Ao = ^^^kBT f‘^^^Ti[ia2f{eu,PF)], (17) 


where /(e„, pf) = T‘2TTiV for positive and negative Mat- 
subara frequencies, respectively [see Eq. (|^]. A < 0 
is the electron-phonon coupling constant and Cc is the 
high-energy cut-off of the order of the Debye frequency. 
Equation 0 in our case simplifies to 


Aq In = 27rfc_B'r Y 

.„>o 


Dn Aq 


, ( 18 ) 


where T^o is the clean limit superconducting transition 
temperature. In Eq. (181 we have eliminated the cou¬ 


pling constant A in favor of T^n by using the standard 
relatioiP^ 


^ Y _ \ 

-—;- = In- \-2FkBT —. 


(19) 


Besides the magnitude of the order parameter Aq it 
is useful to know the actual superconducting transition 
temperature Tc of the disordered system. Linearizing Eq. 
(181 with respect to Aq we obtain 


G) - G + ssfj; 

Teff = 2Q;^Ugr/dc, dc = {1 + ul — a^Ug)^ + 


( 20 ) 


where Y{x) is the digamma function. Equatio n is 
similar to the famous Abrikosov-Gor’kov formulcpSl, with 
an effective pair-breaking parameter Teff. 

As soon as the self-consistent solution for Aq is found 
we can compute the density of states by using the recipe 
given in Eq. The coherence function becomes. 


70 (e) 


D 


( 21 ) 


while the final result for the density of states reads 


A(e) = 2NFRe 


' l + 7g ' 


= 2AFlm 


E 

y/D^ - E^ 


, ( 22 ) 


where E and D are the real-energy versions of i?„ and 
Dn defined as 


*E'„(e„) —>■ E{e), Dn{en) —t D{e). (23) 

They satisfy 

^ [ul + a^ul + {ul - a^ulf] Ey/D^ - E^ 

^ P + (^0 - a^ulY]Dy/D^ - E^ 

T) = {l + ul- a^ul f{D^ - EY - ia^ulE^. (24) 
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2. Ferromagnetically ordered impurities 


Eq. In this case it reads 


For the case of ferromagnetically ordered impurities, 
see Fig. Sb), there is an exchange field in the system, 
which makes properties of opposite spin quasiparticles 
inequivalent. Thus it is convenient to parametrize the 
coherence functions and impurity self-energies as 


7 = 


_ f 0 It 

-7t 0 


7 = 


0 


7t 


-7t 0 


E- - i ^ 

^imp - I 0 S 


^inip — 


0 

-A; 0 


(25) 


which would simplify to the randomly oriented case con¬ 
sidered before upon setting the opposite spin-components 
equal. Now we introduce 


Ent = £" + = Cn - i^Fus, (26) 

and solve the t-matrix equation (see Appendix 0, by 
analogy with the previous case. Then we introduce the 
spin-dependent self-energies 


. — F I , F I — F^ 

Ent = Cnt + * ^ j Enl = e„^ -|- i - , 

Dnt = Aq -I- A^, Dnt — Aq -T A^, 


in terms of which the solution to the homogeneous Riccati 
equations is given by 


7T(e«) = * 

7t(en) = * 
7t(en) = - 


Dnt 

Dnj 

+ ^ni 

-7t(e«), 7t(en) = - 


■7T(en)> 


(28) 


and the following symmetries hold 

Enx = Enx, Dnx= Dnx, X = {t; i} • (29) 


Spin-dependent self-energies satisfy the following self- 
consistency equations 


Aq In = 2'KkBT ^ Re 


er.>0 


D 


n't 




n't 


D 


71^ 




Aq 


(31) 


Linearizing this equation with respect to Aq we can also 
find the equation for the critical temperature T^, 




1 . r usja + Pdct) 

2 *27rfcBrc dct 

1 . T us{a + Pdci) 

'^^2 ^2-KkBT^ 


^4 


— 2 I u -|- 


1 


-1 


2 

dct = 1 -I- Uq ~ + 2iaus, del ~ '^ct- 


(32) 


Equation (32) could in principle be rewritten in terms of 
digamma functions of complex argument, but the form 
above was used in the actual calculations. 

The (spin-resolved) density of states is computed via 
analytical continuation to the real energy axis, see Eq. 
@, by analogy with the previous case, and the final ex¬ 
pression reads 


Nt^l{e) = Nplm 


E- 


td 




(33) 


Here E^^i and 0^,1 satisfy 

^ {ul - a^ul)E^ =F - E^ 


+ b- 


— Aq + r 




(ug - a^ul)D^ 
X>v 


, x = {t4}. 


(34) 


V^ = il + ul- a^ul)^Dl-El± 2ausE^, 


^nx r 


{ul - a‘^ul)Enx ± iaus^Dl^ + El^ 


Vr, 


= + x = {t,;}, (30) 


n 


nx 


T^nx ~ E Uq a Ug) ^+ ^nx ^ “StallsEnx, 


where = e =F f^Tug. 

Another observable we are interested in for the case 
of polarized impurities is the magnetization in the sys¬ 
tem. In the framework of our model, the expression for 
magnetization read^^. 


M = Mm, 

M = 2{a + l3)rustiBNF —‘^TrkBTfj.BNp (35) 


where the upper (lower) sign refers to spin-up (spin- 
down) quasiparticles. Next, we find the self-consistency 
equation for the order parameter, which is still given by 


X 


E Im 

e„>0 


Ent 


Enj 

\l^nl + ^nl_ 
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The first term corresponds to the normal state contri¬ 
bution, while the second one is the low-energy correc¬ 
tion due to superconductivity. Note that for unpolarized 
impurities the self-energies are spin-degenerate and the 
magnetization vanishes (the first term vanishes after av¬ 
eraging over impurity directions, see Appendix [A|). 


Finally, if Eq. (31) has more than one solution it is nec¬ 
essary to consider the difference between Gibb s free en- 
ergies in the superconducting and normal state^^^MIlIllIl 
in order to determine the physically relevant one. In our 
case the expression for the free energy difference is given 
by (see Appendix [b| 


= ^qNp In —-h ‘2'kN pkgT \ Re 


e„>0 




n 2 


^n\. 


E, 


'rat 




-^rat + ^nt + ^ni_ 


+ ‘IpNpksT 'y ( Re 


£,.>0 


{En-\ — ^n-\)En^ + (Hrat 7 Ao)Ilnt (^rat ~ f-ni)En\. + (^nt ~ ^o)DnX 


sjDl^ + El 


— 2pNpkBT —Re < In 

e „>0 ^ [ 

-In (l + Uq - a^ug)^ + 


nt 


2 2 2 I 

1 ~h Vjq — Ol Zig ~\~ 


Jd"^, + 

y 




^ , z 2 2 
1 -j- Uq — Oi ZXg — 




^rat+£^rat, 


2fQMg£'n; 

'^rat+£^rat, 


(36) 


The free energy difference for the unpolarized case is 
obtained by equating the opposite-spin self-energies (let 
P = 0). The self-energies then satisfy Eq. (16). 


III. RESULTS 


In this section we use the formulas obtained above to 
investigate properties of superconductors with magnetic 
impurities. Eor each observable we compare results of 
self-consistent numerical calculations for the cases of (i) 
randomly oriented impurity magnetic moments, and (ii) 
ferromagnetically ordered magnetic moments. For the 
ordered case, all numerical results presented below were 
obtained by setting P = 1 and |a| = 0.1 in order to 
illustrate the Physics. These parameters should be con¬ 
sidered as fit parameter when comparing theory with ex¬ 
periment. 

Properties of the system depend very weakly on the 
strength of the scalar part of the impurity potential uq. 
As was demonstrated in Ref. , if the superconducting 
order parameter is isotropic the scalar part Uq enters the 
theory only through the energy of the bound state. This 
can be accounted for by introducing an effective exchange 
scattering amplitude {uo,us} —>■ Mg®. Therefore, in all of 
the results presented below we have taken uq = 0 to 
reduce the parameter space. 


A. Superconducting transition temperature 


We start by computing the superconducting transition 
temperature as a function of impurity concentration. We 
note that there are three ways to do it. First, one can 
solve the linearized gap equation, see Eqs. (20) and (32). 
Second, one can search for the value of temperature for 
which the full self-consistency equation [see Eqs. (18) and 
(31)] for the order parameter has zero solution. Third, 
one can find the temperature at which the free energy dif¬ 
ference Eq. (36) changes sign. All methods give the same 
answer if the superconducting phase transition is second- 
order. When the phase transition is first-order, however, 
the only way to determine the physical transition tem¬ 
perature is by computing the free energy difference. 


1. Randomly oriented impurities 


For the case of randomly oriented impurity spins we ob¬ 
serve suppression of the transition temperature which is 
similar to the classic result by Abrikosov and Gor’ko'^®. 
The difference is that the effective pair-breaking parame¬ 
ter Feff now depends on both the impurity concentration 
and the strength of the impurity potential, see Fig. |^and 
Eq. (20). Rescaling the horizontal axis of the figure with 
respect to Tc/2-KkBTco = dc/Sya^Ug (where 7 « 1.78 is 
the Euler constant), corresponding to the critical den¬ 
sity of impurities at which Tc = 0 , one can see that all 
the curves align, see inset of Eig. The sign of a is 
unimportant for this case. 
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TOTTkeTco 

FIG. 3. (Color online) Main plot: Transition temperature as 
a function of impurity density [see Eq. 01 for randomly 
oriented impurity spins. From right to left ws = 3, 5, 9. Inset: 
The same data with the horizontal axis rescaled with respect 
to Fc, critical density at which Tc = 0 (see text). 


2. Ferromagnetically ordered impurities 

Let us now discuss the case of ferromagnetically or¬ 
dered impurities. In Fig. we plot the critical tem¬ 
perature Tc calculated by the three different methods 
described above. For small impurity concentrations all 
three curves coincide and the superconducting phase 
transition is second order. For high enough concentra¬ 
tions we notice the appearance of the second solution, 
which means that there are two possible values of the 
order parameter in the system at low temperatures, as 
we will see below. The discrepancy between the three 
methods of finding Tc is caused by the fact that the su¬ 
perconducting phase transition becomes first-order. This 
is expected since apart from the impurity scattering we 



r/2TTkBTco 


FIG. 4. (Color online) Critical temperature as a function 
of impurity density for ferromagnetically ordered case and 
o > 0. From right to left we have us = 3, 5, 9. Solid lines are 
solutions to Eq. (321. Dashed lines correspond to 5Q,{T) = 0, 
see Eq. (361. Dotted lines are found from Eq. (311 requiring 
Ao(r) = 0. 



r/2TTkBTco 

FIG. 5. (Color online) Transition temperature as a function 
of impurity density for ferromagnetically ordered case with 
Us = 3. The red (right) curve corresponds to a < 0, while 
the black (left) one is for a > 0, see Fig. 


also have a background Zeeman f ield in our model, see 
Eq. ([^. It is well knowrP^I^SllZIlll ^j^^t for high exchange 
fields the order of the superconducting phase transition 
changes to being first order. For an ordinary Zeeman 
term due to an external in-plan magnetic held, this is 
commonly known as the Clogston-Chandrasekhar limit. 
This means that the order parameter for temperatures 
close to Tc is not small, and it is not allowed to lin¬ 
earize the self-consistency equation (311 to hnd Tc- The 
region in parameter space where this happens is the re¬ 
gion where the solution to the linearized gap equation 
displays a back-bend (see solid lines in Fig. |^. After 
solving the full (non-linear) self-consistency equation, we 
obtain the result depicted by the dotted lines in Fig. 
But, in order for the superconducting phase to exist it 
has to be more energetically favorable than the normal 
one, which can be checked by computing the Free en¬ 
ergy difference. Searching for the temperature at which 
this condition fails, we obtain the physically correct so¬ 
lution depicted by the dashed lines in the hgure. We 
have to mention that in our paper we assume a spatially 
constant and homogeneous order parameter, which cor¬ 
responds to canonical Cooper pairing of electrons with 
equal and opposite momenta. However, there is a possi¬ 
bility of having a solution corresponding to Cooper pairs 
with finite cente r-of-m ass momentum, which is known as 
the FFLO phasd^^lZOJ q^j- Q^gg i^j^jg solution could ex¬ 
ist, in principle, and it would lie in between the dashed 
and dotted lines in Fig. This is, however, out of scope 
of the present paper. 

The results discussed above referred to the case of anti¬ 
ferromagnetic interaction of impurity spins with itinerant 
electrons, a > 0, see Fig. Let us now briefly discuss 
what happens if the local exchange interaction is ferro¬ 
magnetic, a < 0. In Fig. we compare the two cases 
and we can see that they look quite similar. For the pa¬ 
rameters chosen in the figure, the critical temperature 
for q: < 0 is always higher than for a > 0. In order to 
understand why it happens one has to remember that 
























r/2TTkBTco 



T/Tco 

FIG. 6. (Color online) Order parameter as a function of im¬ 
purity density (a) and temperature (b) for randomly oriented 
impurity spins, (a): From right to left Us = 3,5,9 and T = 
0.01 TcO- (b): From top to bottom V/2'KkBTco = 0.1,0.2,0.3 
and Ms = 3. 



r/2TTkBTco 


FIG. 7. (Color online) Order parameter (a) and the free en¬ 
ergy difference (b) as a function of impurity density for fer- 
romagnetically aligned impurity spins and a > 0. From right 
to left Ms = 3,5,9 and T = 0.01 Tco. Two different solutions 
to Eq. (311 and the corresponding free energy difference are 
plotted with solid and dotted lines. Thin vertical lines mark 
the density at which the free energy difference crosses zero. 


the suppression of the order parameter (and, as a con¬ 
sequence, of the critical temperature) as a function of 
impu rity d ensity is caused by the growing band of YSR 
stated^oissHU] shrinking the energy gap in the spectrum. 
Discussing the density of states below we will show that 
changing the sign of a changes the spin-polarization of 
the YSR states to the opposite one. Since we also have a 
Zeeman-like shift in our model, see Eq. ([^, which does 
not depend on the value of a, the decrease of Tc depends 
on how soon the impurity band meets the quasiparticle 
continuunPS. Therefore, one can find a set of parame¬ 
ters when, in contrast to the results in Fig. the critical 
temperature for a < 0 is smaller than for a > 0. 


the reduction of Aq as a function of impurity concen¬ 
tration is slower than the corresponding reduction of the 
critical temperature, see Fig. Therefore if one plots 
Ao/ksTc as a function of density it seems that the order 
parameter grow^^(not shown in the figure). In Fig. |^b) 
we plot the order parameter as a function of temperature 
for three different impurity densities. All the plots have 
similar shape and simply reflect the gradual reduction of 
Tc and Aq. 


2. Ferromagnetically ordered impurities 


B. Order parameter 

1. Randomly oriented impurities 


For randomly oriented impurity spins the order param¬ 
eter behaves similarly to the Born limit considered by 
Abrikosov and Gor ’ko^^iSl. 

As can be seen in Fig. |^a). 


Let us now consider the case of ferromagnetically or¬ 
dered impurity spins. In Fig. 0^) we plot the order 
parameter as a function of impurity density for different 
values of us- For low densities there is only one solution 
to Eq. (31), which decreases monotonically. However, 
for large enough densities the second solution emerges, 
depicted by the dotted lines in Fig. [^a). In order to de¬ 
termine the physically relevant one we plot the difference 
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FIG. 8. Order parameter (a)-(c) and Gibbs free energy difference (d)-(f) as a function of temperature for ferromagnetically 
aligned impurity spins with us = 3 and r/2'KkBTco equal to 0.04 for (a) and (d), 0.045 for (b) and (e), 0.055 for (c) and (f). 
The vertical line in (f) depicts the temperature at which the free energy difference crosses zero. The inset in (d) shows a zoom 
on the dotted line in the main plot. The dotted lines in (b)-(c) depict the second solution to the order parameter equation | |31[ ) 
and the corresponding free energy difference in (e)-(f). 


of free energies in superconducting and normal states, 
see Fig. As can be seen from the figure, the new 

solution has a non-negative free energy difference, which 
means that it is not energetically favorable and thus does 
not realize in practice. A similar conclusion was found 
in Ref. [48) . where the superconductor to normal metal 
phase transition in a pure Zeeman field was analyzed. 
Another interesting point is that for even higher densi¬ 
ties both solutions become energetically unfavorable and 
the system is no longer superconducting (in fact, it be¬ 
comes a ferromagnet due to impurity ferromagnetisirP^ . 
The transition point at which it happens is depicted by 
thin vertical lines in Fig. [^b). If we change the sign 
of a, corresponding to a different type of the local ex¬ 
change interaction with impurities (see Fig. [2 1, a similar 
behavior is observed, thus we do not show it here. 

To demonstrate the temperature dependence of the or¬ 
der parameter, we choose three different values of F. In 
Figs, [^a)-(c) we plot the order parameter as a function 
of temperature for us = 3 and F = 0.04, 0.045 and 0.055, 
respectively (see black lines in Figs. |^and[^. In Figs. 
[^d)-(f) we show the corresponding free energy differ¬ 
ence for each order parameter solution (in case there are 
more than one). All plots demonstrate that there is only 
one physically relevant solution (with the biggest value of 
Aq), for which Sfl < 0. Moreover if the impurity density 
is far enough from the point where 517 = 0 [thin vertical 
line in Fig. [^b)], the order parameter goes gradually to 
zero indicating that the phase transition is of second or¬ 
der, see Figs. |^a)-(b). Otherwise, the order parameter 
goes abruptly to zero at the critical point where the free 
energy difference becomes positive and the phase transi¬ 


tion is hrst-order, as in Figj^c). Changing the type of 
local exchange scattering off impurities by inverting the 
sign of a does not alter qualitatively the results discussed 
above and it is not shown here. 


C. Density of states 


1. Randomly oriented impurities 


For the case of unpolarized impurities opposite-spin 
densities of states are the same and we therefore plot 
the total density of states as a function of energy in Fig. 
1^ Figure §a) demonstrates evolution of the density of 
states with varying impurity potential strength us, while 
Fig. I^b) shows its evolution with impurity density F. 
As one can see, the impurity potential strength sets the 
position of the YSR bands inside the gap, keeping the 
total number of YSR states fixed. On the other hand 
the density of impurities sets the size of the bands, by 
increasing the overall number of YSR states. Figure|^b) 
also shows that randomly oriented impurity spins play a 
pair-breaking role in the system, decreasing the electron- 
hole coherence. The latter result was also obtained by 
Abrikosov and Gor’ko^ff^in the Born limit. They demon¬ 
strated that the order parameter and the energy gap in 
the spectrum are, generally speaking, two different quan¬ 
tities and there is a range of impurity densities where the 
so-called gapless superconductivity emergeJi^. 
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£/Ao(r=o) 


FIG. 9. (Color online) Density of states for randomly ori¬ 
ented impurity spins, (a): Evolution with us. Black, red and 
blue lines correspond to us = 2,4,8 and V/2'KkBTco = 0.01. 
(b) Evolution with F. Black, red and blue lines correspond 
to T/2'KkBTcQ = 0.05,0.1,0.2 and us = 5. For both plots 
temperature is T = O.OlTco. 


2. Ferromagnetically ordered impurities 


For the case of ferromagnetically ordered impurities, 
opposite-spin densities of states are distinct and we plot 
them separately in Fig |10| The most important differ¬ 
ence from the previous case is that YSR impurity bands 
are now spin-polarized, which allows for a non-zero mag¬ 
netization in the system, as we will see below. Another 
complication arises from the background homogeneous 
exchange field generated by the impurities [see Eq. ([^], 
which shifts the opposite-spin spectra with respect to 
each other. The latter circumstance is important for un¬ 
derstanding the weakening of superconductivity in this 
case. Indeed, if one increases impurity strength ug, the 
energy shift due to the exchange field drives the system 
to the normal s tate s imilarly to the effect of a pure Zee- 
man interactioiP^I^. On the other hand, increasing the 
density of impurities results in a simultaneous growth of 
the impurity band and a shift due to the exchange field, 


which leads to faster decrease of superconducting prop¬ 
erties than in the unpolarized case. Comparing Fig. ©a) 
and Fig. one can see that allowed densities of mag¬ 
netic impurities for this case are one order of magnitude 
smaller than those for the unpolarized case. 

Next, we compare how the type of local exchange 
scattering off impurities (see Fig. changes the spec¬ 
tral properties of the system. For the case of anti¬ 
ferromagnetic scattering (black lines in Fig. 101, a > 0, 
the spin-up YSR impurity subband splits off from the 
quasiparicle continuum at negative energies (solid lines), 
while spin-down sets in at positive energies (dotted lines). 
At the same time for ferromagnetic interaction (red lines 
in Fig. 10), a < 0, the opposite-spin YSR subbands ex¬ 
change places. As can be seen from Figs. [To|[a)-(c), the 
Zeeman-like shift due to the background exchange field 
counteracts the impurity band shift imposed by varying 
us for a < 0, while the two effects cooperate for a > 0. 
In Figs. -(f) we observe that the background ex¬ 

change field shifts the growing impurity subbands in the 
opposite directions for a > 0 and a < 0. All these fea¬ 
ture are responsible for the differences in the transition 
temperature Tc in Fig. 


D. Magnetization and Sakurai phase transition 


When spins of magnetic impurities are ferromagneti¬ 
cally aligned, the single-particle spectrum of the system 
is not spin-degenerate, as we have seen in the previous 
section. This manifests itself as an imbalance between 
occupations of opposite-spin subbands, resulting in net 
magnetization. Moreover, since the in-gap YSR impurity 
bands are spin-polarized (see Fig. 10), we demonstrate 


below how one can observe the signature of the quantum 
phase transition first discussed by Sakurai in Ref. [57]. In 
this seminal paper it was shown that in presence of a sin¬ 
gle YSR magnetic impurity inside a superconductor there 
is a quantum phase transition in the ground state of the 
system when the impurity-induced YSR state crosses the 
Fermi level. By increasing the effective coupling strength 
C = nNpJS, where J is the exchange interaction con¬ 
stant and S is the impurity spin, the YSR state moves 
from one side of single-particle continuum to the opposite 
one, crossing zero when C = 1- The qualitative physical 
picture of the phase transition is as follows. Since the 
YSR state is spin-polarized, when it gets occupied a sin¬ 
gle itinerant electron from the superconductor is bound 
to the impurity site in a singlet or triplet state depending 
on the nature of the local exchange interaction, see Fig. 
m At the same time its time-reversed mate is left with an 
uncompensated spin and the ground state of the syste m 
in this case always contains a single quasiparticl^^HHl, 

In terms of the parameters of our model, we have 
irNpJS = aus and, in addition, there is a finite den¬ 
sity of impurities n = irNpT, instead of a single impu¬ 
rity. The parameter a has a meaning of the tunneling 
amplitude onto the impurity site, which means that the 
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(a) 



(d) 




£/Ao(r=o) £/Ao(r=o) £/Ao(r=o) 


FIG. 10. (Color online) Spin-up (full line) and spin-down (dotted line) densities of states for ferromagnetically aligned impurity 
spins. Black lines correspond to a > 0 and red ones to a < 0. (a)-(c) show the effect of increasing us ~ 2,4,8 with 

r 12'KkBTc.o = 0.01. (d)-(f) show the effect of increasing F j^'KksTca = 0.005,0.01, 0.02 with us = 5. Temperature isT = O.OlTco. 


parameters that can be controlled, in principle, are ug 
and r. Using Eq. (35) we compute the magnetization 
M. By definition, the total magnetic moment in the sys¬ 
tem is = MA, where A is the volume of the system. 
On the other hand, the total magnetic moment is re¬ 
lated to the total spin of the system (in units of h) S via 
A4 = —g^sS, where g is the quasiparticle g-factor and 
/is is the Bohr magneton. Then we can introduce the 
average spin per magnetic impurity s = S/N via 


M = —gTrNpTg.BS. (37) 


In Fig. [TT|(a) we plot s as a function of impurity 
strength us- As we know from the discussion of the den¬ 
sity of states above, varying us changes the position of 
the YSR impurity band inside the gap, see Figs. [^a)- 
(c). Therefore by tuning this parameter it is possible to 
make the impurity band cross the Fermi energy, which 
would mean that each magnetic impurity binds a single 
electron to itself, according to the qualitative picture of 
the phase transition discussed above. Then the average 
spin per impurity which is left uncompensated after the 
phase transition should be s = =f 1/2 for sgn(a) = ±1. 
As can be seen from Fig. [TT|(a), there is indeed a jump in 
s exactly equal to this value, but with a finite slope de¬ 
termined by the impurity density T. In addition, there is 


a background magnetization which makes the magnitude 
of s slowly increasing before and after the transition. It is 
determined by the difference in the densities of states in 
the (negative-energy) quasiparticle continuum, while the 
jump is due to the YSR impurity band. This difference is 


hardly seen in Fig. 10 but the more the YSR band splits 
off, the more it deforms the continuum it originated in, 
and the bigger the difference between the opposite-spin 
densities of states becomes. This additional feature was 
left out in the qualitative picture proposed by SakuraP^. 
The free energy difference in Fig. ©b) shows that the 
system remains superconducting after the phase transi¬ 
tion. If a < 0, the results look similar, except that s > 0. 

In Fig. [^we perform a similar analysis, but we vary 
impurity density instead of their strength. This case is 
more suitable for experiments since controlling the den¬ 
sity is easier. Moreover, in this case the phase transition 
is much more clearly observed. Indeed, in Fig. [^a) we 
plot the difference between the average spin per impurity 
in superconducting s and normal sn states. The latter 
is given by the first term on the rhs of Eq. (35). One 


can see a clear jump in the figure with amplitude equal 
to —1/2, as expected. We already know that varying T 
makes the YSR impurity band grow and shift (due to the 
background exchange field) at the same time, see Figs, 
©d) -(f). Thus, there is a limited range of parameters 














































12 



0 2 4 6 8 10 12 14 16 


Us 


FIG. 11. (a) Average spin per impurity s [see Eq. (37l] as a 
function of us for a > 0. The impurity density F /2nkBTcQ is 
equal to 0.001 (solid line), 0.005 (dashed line) and 0.01 (dotted 
line), (b) The corresponding Gibbs free energy difference. 


where one can see both plateaus, as in Fig. [T^a). More¬ 
over, because of the pair-breaking effect of the impurities, 
some of the results are not observable because of tran¬ 
sition to the normal state, indicated by the free energy 
difference in Fig. dltb). Finally, if a < 0, in order to see 
the jump associated to the phase transition by varying 
F, one has to start in a state when the transition has oc- 
cured (by choosing high enough us, see Fig. 10). Then, 
gradually increasing F the impurity band gets empty in¬ 
stead of becoming occupied (as for a > 0), and the initial 
plateau at 1/2 evolves to the final one at zero. 


IV. DISCUSSION AND CONCLUSIONS 

Let us now consider in more detail our theoretical 
model and its range of validity. We start from a dis¬ 
cussion of the necessity to take into account the back¬ 
ground magnetic field collectively created by the distri¬ 
bution of magnetic impurities, see Eq. ([^. This ex¬ 
change field makes the interaction of itinerant electrons 
with the magnetic impurities non-local. This is expected 
to hold if the impurity magnetic moment increases, since 
the magnetic field created by a magnetic dipole decays 



r/2TTkBTco 

FIG. 12. (Golor online) (a) Difference in the average spin (per 
impurity) in the superconducting s and normal sjv phases 
as a function of impurity density for q > 0. The impurity 
strength, from right to left, is us = 6,7, 8,9. (b) The cor¬ 
responding Gibbs free energy difference. Thin vertical lines 
depict the critical density at which it crosses zero. 


as 1/r^ as a function of distance r from the dipole cen¬ 
ter. In our model increase of the magnetic moment of 
impurities corresponds to increasing its. An interest- 

^ feature of the t-matrix approximation, see Appendix 
is that the single-impurity t-matrix timp(en) becomes 
spin-independent if its oo (unitary scattering limit), 
see Eq.(30). This is a formal result and it means that 
in this limit, the magnetic impurities would behave just 
like scalar scatterers, which are known to satisfy the An¬ 
derson’s theoreirP^. This would lead to a paradox that 
for the case of ferromagnetically ordered impurities, if 
we increase their strength its, we recover the clean-limit 
results for the order parameter and other properties of 
the systerrP^. The paradox is solved when we include 
the magnetic field generated by the ordered impurity 
moments, see Eq. (§. In this case the order parameter 
gets suppressed long before we achieve the unitary limit. 
Note that this complication is absent when impurities are 
treated in the Born limilPSll 

because then the only free 
parameter is the impurity density (the impurity strength 
is assumed small). 

We would also like to comment on the description of 
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the impurity subsystem. In this work we considered the 
two limiting cases: (i) unpolarized and (ii) ferromagneti- 
cally ordered impurities. We did not include any theoret¬ 
ical description of how the ferromagnetic ordering takes 
place. However, as soon as the transition to the polarized 
case has occured, our results should be valid. A rigorous 
way of solving this problem would require to include, for 
example, a Heisenberg model for magnetic impurities and 
describe the dynamics of the coupled superconductor- 
impurities system. This would allow to consider a tran¬ 
sient regime when the impurities are partially polarized. 
At the same time one would be able to introduce a new 
energy scale /csTcurie, rel ated to the Curie temperature 
of the impurity subsystenP^EIl^ Then, our results for fer- 
romagnetically aligned impurity spins discussed in this 
work are valid for T Tcurie- Another point related to 
the description of impurities which was omitted here is 
the possibility to form clusters. To our knowledge, this ef¬ 
fect cannot be automatically included into the quasiclas- 
sical theory described in this paper, and it would require 
a separate treatment. Allowing the magnetic impurities 
to cluster would smear out the sharp edges of the YSR 
impurity bands via the so-called Lifshitz tail^^. 

Before conclusion we would like to comment on the or¬ 
der of magnitude of the impurity magnetic moment used 
in this work, corresponding to the value of us- Assuming 
a quadratic energy dispersion of electrons in the normal 
state and taking their effective mass equal to the bare 
electron mass, we can estimate that for a superconduc¬ 
tor with Fermi energy Ep — W eV, the typical size of 
an impurity is ^ 1 — 8 nm, and its magnetic moment 
~ 8.5 X 10^ (for us = !)• This coincides with the 
typical size of s mall ferromagnetic islands used in mod¬ 
ern experiment^i^^^. On the other hand it justifies our 
treatment of the impurity spins as classicaP^. 

In conclusion, we have studied thermodynamic proper¬ 
ties of a superconductor with a finite density of magnetic 
impurities, described within a generalized Yu-Shiba- 
Rusinov model (self-consistent t-matrix approximation). 
When the impurity spins are randomly oriented, most of 
the results are similar to the ones obtained by Abrikosov 
and Gor’ko's^ within the first-order Born approximation. 
The only difference is that the YSR impurity bands are 
split off from the quasiparticle continuum. For the case 
of ferromagnetically ordered impurities we argue that it 
is necessary to include a background magnetic exchange 
field created by their spins. We have found that in this 
case the superconducting transition changes from second 
order to first order as the impurity strength is increased. 
At a critical impurity strength, superconductivity disap¬ 
pears. We have shown that the signature of the quantum 
phase transition of the system ground state due to the 
YSR impurity band crossing the Fermi energy is a drop 
in the magnetization. We emphasize that the initial idea 
of this phenomenon was put forward by SakuraP^ for a 
single impurity, while we have shown how it manifests 
itself for a finite impurity density. 
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Appendix A: Impurity self-energy in t-matrix 
approximation 

In terms of the Matsubara Green’s function, the impu¬ 
rity self-energy in t-matrix approximation is given by^ 

—matr (^n 5 Pi^) — ^^imp(Cn; Pi^ 1 Pi^)5 ('^^) 

where the single-impurity t-matrix timp satisfies 

iimp{en, Pf, Pf) = ^(Pf, Pf) (A2) 

f dOp// 

'^(PF, PF)5(£ni PF)^imp(en, Pfj Pf)- 


Here 'C(pF)P^) is the matrix element of the impurity 
potential between the quasiparticle states with momenta 
Pi? and Pp on the Fermi surface (computed in the nor¬ 
mal state of the system). Equation (A2) is usually rep¬ 


resented diagrammatically as a sum of diagrams for ar¬ 
bitrary number of quasiparticle scatterings on a given 
impurity, see Fig. |13| In this paper we consider only 
s-wave scattering off impurities, i.e. u(pfiPf) inde¬ 
pendent of momenta. For the two models of magnetic 
impurities described in the main text, see Eqs. (§-§, 
the matrix element v is 


V 0 
0 V* 


V = vq + ausm • cr. 


(A3) 


1. Self-energy for randomly oriented impurities 


It is worth mentioning that for deriving Eq. (A21 one 


has to perform averaging of the Dyson equation for the 
(full microscopic) propagator over the impurity positions 
defined as [see Eq. 0 ] 


imp. pos. 


N 

n 

1=1' 


A 


(A4) 


X X x, 

timp(en, PFj Pf) = + “•-^ ^ 

X 


w(Pf,Pf)= 


g(e„,pF) = 


FIG. 13. Diagrammatic representation of equation (|A2l). 
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where integration is performed over the system volume 
A. For the case of unpolarized magnetic impurities, be¬ 
sides averagi ng o ver impurity positions, one also has to 
average Eq. (A2) over the magnetic moment directions. 


This is performed by parameterizing the unit vector m 
in spherical coordinates and defining 


*)spin dir. — 


dfir^ 


47r 


(A5) 


Appendix B: Free energy functional 

In this section we briefly describe how to compute the 
difference between the Gibbs free energies in supercon¬ 
ducting and normal states. The free energy is a func¬ 
tional of the quasiclassical propagator and self-energies. 
So, we define 


snig, h, T] = r!s[g, h, T] - r!^[ 5 , h, T], (Bl) 


Thus, using equations and ( |A2[ ) one can obtain 

the self-energy 


which means that for 6fl < 0 , the superconducting state 
i s m ore energetical ly favor able. In order to derive Eq. 
(36) we follow Refsand write down 


^imp(Cn) — ^(^imp (Cn))spin dir. 7 


with matrix elements written in Eq. (13). 


(A6) 


2. Self-energy for ferromagnetically ordered 
impurities 


For the case of ferromagnetically ordered magnetic im¬ 
purities we can choose the coordinate system in spin 
space such as nij = m = (0,0,1). Since in this case, 
apart from the local scattering by the impurities, we also 
have a background magnetic field in the system [see Eq. 
([^], the impurity self-energy consists of two parts 

^imp (e„) = /dnvsCTsi + n£imp(e„), (A7) 


where the first term has a form of Zeeman interaction, 
while the second one is obtained by solving Eq. (A2). 


ns[g,KT]=nN[9,KT = 0] 


-In, 


drir- 


Att 


-ksT ^ Tr h(pF,e„) 5 (k, e„ 

kriKCc 


d^ln 


-Gq (k,e„)-b h(pF,e„) -b (54> [ 5 ]. (B2) 


Here e„) = ie„f 3 —where is a single-particle 

spectrum in the normal state (calculated with respect to 
the Fermi energy Ep). [g] = [g] — [ 5 ], where 

$ [ol is a functional which generates the perturbation ex- 
pansion for the skeleton self-energy diagram d”*^ * ”” ! We 
note that the log-term on the rhs of Eq. (B2) contains a 
finite temperature contribution to the normal state free 
energy HAr[g, h, T > 0 ], which has to be subtracted when 
com putin g the integral^. Equation (36) is obtained from 
Eq. (B2) using the quasiclassical self-energy h appropri¬ 
ate to our model and assuming quadratic energy spec¬ 
trum in the normal state ^k = h^k^/2m* — Ep, where 
m* is the effective mass. 
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